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FOREWORD 


Acid  deposition  occurs  when  acidifying  pollutants  emitted  from  anthropogenic  and  other 
processes  undergo  chemical  reactions  in  the  atmosphere  and  fall  to  the  earth  as  wet 
deposition  (rain,  snow,  cloud,  fog)  or  dry  deposition  (dry  particles,  gas).  Acidic 
pollutants  can  be  transported  long  distances  in  the  atmosphere  from  their  sources  and 
eventually  be  deposited  in  ecosystems  over  broad  regional  scales  and  in  locations  far 
from  the  emission  sources. 

Dry  deposition  is  generally  more  a  local  problem  than  wet  deposition.  Direct 
measurement  of  dry  deposition  rates  is  difficult.  Dry  deposition  depends  on  many 
factors,  including:  meteorological  conditions,  characteristics  of  the  pollutants  being 
deposited  (e.g.  different  gaseous  chemical  and  particle  size),  and  characteristics  of  the 
surface  on  which  deposition  occurs. 

The  most  accepted  and  common  method  for  estimating  dry  deposition  is  the  so-called 
"inference  method."  The  inferential  method  is  a  combination  of  measurement  and 
modeling  that  involves  indirect  estimation  of  dry  deposition  rates  on  the  basis  of 
routinely  measured  air  concentrations  and  meteorological  parameters.  The  method  is 
based  on  an  assumed  steady-state  relationship  F  =  Vd  C,  where  the  dry  deposition  flux  or 
rate  (F)  is  a  product  of  the  dry  deposition  velocity  (Vd)  and  the  concentration  (C)  of  an 
airborne  pollutant. 

A  series  of  studies  have  been  initiated  by  AENV  to  evaluate  the  inference  method  and 
search  for  the  most  suitable  and  simple  model  for  deposition  rate  estimations  in  Alberta. 
This  report  documents  the  second  study  in  the  series.  Titles  for  the  reports  of  the  other 
studies  are:  ''Review  and  Assessment  of  Methods  for  Monitoring  and  Estimating  Dry- 
Deposition  in  Alberta",  and  "Refinement  Study  of  Dry  Deposition  Inference  Method 
Used  in  Alberta    It  is  anticipated  that  once  all  necessary  information  is  gathered,  an 
Alberta  protocol  for  dry  deposition  measurement  will  be  prepared. 


Lawrence  Cheng,  Ph.  D. 
Air  Policy, 

Climate  Change,  Air  and  Land  Policy  Branch 
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SUMMARY 


Currently  there  is  no  standard  method  for  field  measurement  and  estimation  of  dry 
deposition  of  acidifying  pollutants  in  Alberta.  Alberta  Environment  is  pursuing 
development  of  an  inexpensive  standard  method.  This  may  involve  a  combination  of 
continuous,  intermittent,  and/or  passive  field  measurements  of  acidifying  pollutants  and 
using  inference  methods  for  estimating  sulphur  and  nitrogen  species  dry  deposition.  An 
overall  objective  of  this  study  was  to  continue  evaluating  inferential  methods  for 
estimating  dry  deposition  rates. 

Dry  Deposition  Methods  Evaluation 

Three  different  inference  methods  used  in  Canada  were  evaluated  for  estimating  sulphur 
and  nitrogen  species  dry  deposition.  These  methods  were  developed  by  Alberta 
Environment  (AENV),  Environment  Canada  (ENVC),  and  the  Wood  Buffalo 
Environmental  Association  (WBEA).  A  specific  objective  was  to  use  these  inference 
methods  to  compare  contributions  of  gaseous  SO2,  NO2,  HNO3,  HNO2,  and  particulate 
forms  of  sulphur  (S04^')  and  nitrogen  (NO3'  and  NH/)  in  total  sulphur  and  nitrogen 
deposition  and  their  potential  relationships  among  each  method.  Another  specific 
objective  was  to  identify  and  recommend  a  preferred  method  (or  methods)  that  could  be 
routinely  used  for  estimating  sulphur  and  nitrogen  species  dry  deposition  rates.  The 
following  findings  are  noted: 

1.  Complete  sulphur  and  nitrogen  gaseous  and  particulate  species  and 
meteorological  data  for  a  one-year  period  (2003)  from  the  Wood  Buffalo 
Environmental  Association  Fort  McKay  air  monitoring  station  were  evaluated 
using  the  three  methods.  All  three  methods  showed  a  similar  ability  to  represent 
variation  in  dry  deposition  rates  for  individual  species  from  month  to  month. 

2.  The  WBEA  method  yielded  higher  monthly  deposition  rates  for  most  species 
modeled  (S04^",  NO2,  HNO3,  HNO2,  and  NH4^)  compared  to  the  AENV  and 
ENVC  methods.  Monthly  potential  acid  input  (PAI)  associated  with  dry  species 
for  each  method  is  shown  in  Figure  ES-1.  The  WBEA  method  also  yielded 
higher  PAI  associated  with  dry  species  compared  to  the  AENV  and  ENVC 
methods  during  winter  and  spring.  The  WBEA  method  gave  similar  PAI  values 
associated  with  dry  species  during  summer  and  fall.  A  simpler  computational 
approach  for  aerodynamic  and  surface  resistance  terms  in  the  WBEA  method  may 
partially  explain  differences  observed  during  winter  and  spring.  Further  testing 
would  be  required  to  better  understand  whether  this  is  the  case. 

3.  Annual  total  PAI  associated  with  dry  species  for  2003  was  0. 14  kg  H^/ha/yr 
(AENV  method),  0.15  kg  H^/ha/yr  (ENVC  method),  and  0.20  kg  H"'/ha/yr 
(WBEA  method).  The  higher  annual  rate  observed  for  the  WBEA  method 
relative  to  other  methods  is  consistent  with  higher  monthly  deposition  rates  for 
individual  species.  While  the  AENV  and  ENVC  methods  yielded  comparable 
monthly  PAI  values  associated  with  dry  species,  it  is  observed  that  ENVC  input 
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parameters  have  more  boundary  conditions  (i.e.  cut-off  points  and  lower/upper 
limits)  associated  with  their  use  compared  to  AENV  input  parameters. 


0.04 


Figure  ES-1.  Monthly  potential  acid  input  (PAI)  values  associated  with  dry  species 
for  the  period  January  to  December  2003  at  WBEA  Fort  McKay  air 
monitoring  station  AMS  #1  (expressed  in  units  of  kg  H+/ha/month). 

(ENVC  =  Environment  Canada  method;  AENV  =  Alberta  Environment  method; 
WBEA  =  Wood  Buffalo  Environmental  Association  method) 


4.  The  ENVC  method  is  more  complicated  given  that  it  has  the  greatest  number  of 
input  parameter  requirements.  A  number  of  additional  assumptions  for  some  of 
the  input  parameters  are  also  required  before  modeling  can  be  performed.  The 
AENV  method  provided  comparative  results  to  the  ENVC  method.  The  AENV 
method  is  less  complicated  compared  to  the  ENVC  method,  and  it  requires  less 
input  data  and  computing  time.  The  WBEA  method  is  simple  and  straight- 
forward to  use.  However,  this  method  yielded  higher  monthly  deposition  rates  for 
most  species  modeled  compared  to  the  other  methods  using  the  2003  dataset. 

5.  The  dataset  used  to  evaluate  the  methods  was  short  -  one  year  -  and  month  to 
month  variation  or  lack  of  variation  in  species  deposition  rates  observed  using 
each  method  does  not  necessarily  provide  an  indication  of  what  may  happen  in 
other  years.  Although  only  a  one-year  dataset  was  available  for  evaluation. 
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results  support  that  the  AENV  method  offers  a  practical  choice  for  routine 
computation  of  sulphur  and  nitrogen  species  dry  deposition  in  Alberta.  Further 
testing  using  different  datasets  elsewhere  in  Alberta  is  recommended  to  support 
this.  If  the  ENVC  method  can  be  shown  to  provide  comparable  results  to  the 
AENV  method  using  different  datasets,  it  would  offer  a  logical  choice  for 
periodic  use  as  a  check  on  the  AENV  method. 

6.  Evaluation  of  contributions  of  SO2  and  NO2  deposition  in  total  sulphur  and 
nitrogen  species  deposition  at  the  Fort  McKay  site  indicated  that  about  82%  is  in 
the  form  of  gaseous  SO2  with  the  remainder  as  particulate  804""  (18%).  These 
results  are  comparable  to  annual  deposition  loadings  of  sulphur  species  during 
1998  to  2002  at  the  former  Alberta  Environment  Beaverlodge  acid  deposition 
monitoring  site  (80%  in  the  form  of  gaseous  SO2  and  20%  as  particulate  S04^'). 

Estimates  of  annual  nitrogen  species  deposition  at  the  Fort  McKay  site  indicated 
that  -51%  was  from  NO2  with  the  remainder  as  HNO3  and  HNO2  (35%)  and 
particulate  NH4^  and  NO3'  (<14%).  Estimates  of  annual  nitrogen  species 
deposition  at  the  Beaverlodge  site  during  1998  to  2002  showed  that  about  35  to 
50%  was  from  NO  and  NO2  with  the  remainder  as  HNO3  and  HNO2  (40  to  60%) 
and  particulate  NH4^  and  NO3"  (<10%). 

7.  Future  work  evaluating  these  inference  methods  should  consider  the  influence  of 
assigning  additional  boundary  conditions  to  AENV  model  input  parameters  on 
variability  associated  with  model  outputs.  Specifically,  assigning  boundary 
conditions  to  AENV  input  parameters  consistent  with  default  ENVC  model 
assumptions  should  be  investigated  to  better  understand  whether  these  changes 
have  the  ability  to  improve  AENV  model  performance  relative  to  the  ENVC 
model. 
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1.0  INTRODUCTION 


Acid  deposition  occurs  when  acidifying  pollutants  emitted  from  anthropogenic  and  other 
processes  undergo  complex  chemical  reactions  in  the  atmosphere  and  fall  to  the  earth  as 
wet  deposition  (rain,  snow,  cloud,  fog)  or  dry  deposition  (dry  particles,  gas).  The  main 
chemical  precursors  leading  to  acidic  pollutants  are  sulphur  dioxide  (SO2)  and  oxides  of 
nitrogen  (NOx).  Reactions  of  these  pollutants  with  water,  oxygen,  carbon  dioxide,  and 
sunlight  in  the  atmosphere  produce  acidic  pollutants,  e.g.  sulphuric  acid  (H2SO4),  nitric 
acid  (HNO3),  and  nitrous  acid  (HNO2).  These  and  other  acidic  pollutants  can  be 
transported  long  distances  in  the  atmosphere  from  their  sources  and  eventually  be 
deposited  in  ecosystems  over  broad  regional  scales  and  in  locations  far  from  the  emission 
sources. 

Dry  deposition  refers  to  removal  of  aerosol  pollutants  through  eddy  diffusion  and 
impaction,  large  particles  through  gravitational  settling,  and  gaseous  pollutants  through 
direct  transfer  from  air  to  water  via  gas  exchange.  Dry  deposition  involves  acidic  sulphur 
and  nitrogen  pollutants  (gases  or  particles)  from  the  atmosphere  being  retained  by  the 
earth's  surface.  At  the  same  time,  co-deposition  of  base  cations  (e.g.  Na^,  Mg^"^,  Ca^^  and 
K^)  results  in  a  reduction  of  the  amount  of  deposited  acidity. 

Potential  acid  input  (PAI)  provides  a  method  of  representing  total  acidic  deposition 
(Cheng  et  al.,  2001).  PAI  includes  both  wet  and  dry  deposition.  PAI  is  calculated  by 
subtracting  neutralizing  capacity  (base  cation  deposition)  from  deposition  of  acidic 
substances  (e.g.  sulphur  and  nitrogen  species).  Cheng  et  al.  (2001,  1997)  provide  a 
detailed  description  of  the  estimation  of  total  PAI.  The  PAI  method  does  not  include 
processes  that  remove  acidity  from  the  earth's  surface  (leaching,  runoff,  etc.). 

1.1      Objectives  of  Study 

Currently  there  is  no  standard  method  for  the  field  measurement  and  estimation  of  dry 
deposition  of  acidifying  pollutants  in  Alberta.  Alberta  Environment  is  pursuing 
development  of  an  inexpensive  technique  for  measuring  and  using  inference  methods  for 
estimating  sulphur  and  nitrogen  species  dry  deposition.  WBK  (2008)  reviewed  current 
approaches  used  for  measuring  and  estimating  dry  deposition. 

An  overall  objective  of  this  current  study  was  to  continue  evaluating  inferential  methods 
for  estimating  dry  deposition  after  WBK  (2008).  The  purpose  was  to  further  investigate 
inferential  methods  that  can  be  used  by  government  and  industry  in  Alberta  for  improved 
acid  deposition  management.  Three  different  inference  methods  used  in  Canada  for 
estimating  sulphur  and  nitrogen  species  dry  deposition  were  evaluated.  A  specific 
objective  was  to  compare  contributions  of  gaseous  SO2,  NO2,  HNO3,  HNO2,  and 
particulate  forms  of  sulphur  and  nitrogen  in  total  sulphur  and  nitrogen  dry  deposition  and 
identify  potential  relationships  of  each  method.  Another  specific  objective  was  to 
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identify  and  recommend  a  preferred  inference  method  (or  methods)  that  could  be 
routinely  used  for  estimating  sulphur  and  nitrogen  species  dry  deposition  rates  in  Alberta. 
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2.0    DRY  DEPOSITION  METHODS  EVALUATION 


The  Wood  Buffalo  Environmental  Association  (WBEA)  collected  complete  sulphur  and 
nitrogen  gaseous  and  particulate  species  and  meteorological  data  at  the  Fort  McKay  air 
monitoring  station  -  AMS  #1  -  during  2003.  These  data  were  used  to  evaluate 
relationships  of  sulphur  and  nitrogen  species  deposition  using  three  inference  methods 
applied  in  Canada.  These  methods  were  developed  by: 

•  Alberta  Environment  (Cheng  et  al.,  2001). 

•  Environment  Canada  (Zhang  et  al.,  2003a,  2003b,  2002a,  2002b,  2001a,  and 
2001b). 

•  WBEA  Terrestrial  Environmental  Effects  Monitoring  (TEEM)  Committee 
(EPCM,  2002,  2000). 

A  specific  objective  of  the  evaluation  was  to  compare  contributions  of  gaseous  SO2,  NO2, 
HNO3,  HNO2,  and  particulate  forms  of  sulphur  (S04^")  and  nitrogen  (NOs'  and  NH/)  in 
total  sulphur  and  nitrogen  deposition  and  their  potential  relationships  among  each 
method.  Another  specific  objective  of  the  evaluation  was  to  identify  and  recommend  a 
preferred  inference  method  (or  methods)  that  could  be  routinely  used  for  estimating 
sulphur  and  nitrogen  species  dry  deposition  in  Alberta. 

2.1      Background  on  Inference  Method 

Dry  deposition  is  generally  far  more  a  local  problem  than  wet  deposition.  Estimating  dry 
deposition  rates  is  more  difficult.  Dry  deposition  depends  on  many  factors,  including: 
meteorological  conditions,  characteristics  of  the  pollutants  being  deposited  (e.g.  particle 
size),  and  characteristics  of  the  surface  on  which  deposition  occurs  (US  EPA,  2001). 

A  common  approach  to  indirectly  estimate  dry  deposition  rates  is  on  the  basis  of 
routinely  measured  air  concentrations  and  meteorological  parameters.  Continuous  and/or 
integrated  measurement  techniques  are  used  to  record  the  concentrations  of  atmospheric 
pollutants  and  continuous  measurement  techniques  are  used  to  record  meteorological 
parameters.  These  parameters  are  used  to  estimate  dry  deposition  using  the  most 
common  method  -  the  inference  method  (Cheng  et  al.,  2001;  Brook  et  al,  1999a,  1999b; 
Clarke  et  al.,  1997). 

The  inference  method  involves  indirect  estimation  of  dry  deposition  rates  on  the  basis  of 
routinely  measured  air  concentrations  and  meteorological  parameters.  The  method  is 
based  on  an  assumed  steady-state  relationship  F  =  Vd  C,  where  dry  deposition  flux  or  rate 
(F)  is  a  product  of  the  dry  deposition  velocity  (Vd)  and  the  concentration  (C)  of  an 
airborne  pollutant.  Vd  is  estimated  on  the  basis  of  resistance  models  and  is  defined  as  the 
inverse  of  the  sum  of  multiple  resistance  factors  (aerodynamic  resistance  (Ra),  boundary- 
layer  resistance  (Rb),  and  surface  resistance  (Rc))  (Wesely  and  Hicks,  2000,  1977): 
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(1) 


Figure  1  illustrates  the  relative  locations  where  dry  deposition  resistance  factors  R;,,  Rh. 
and  Rc  apply  near  a  surface. 

Aerodynamic  Resistance  (Ra).  A  shallow  sub-layer  occurs  next  to  the  ground  that  is 
within  the  atmospheric  constant  flux  layer.  The  depth  of  this  layer  is  in  terms  of  meters 
(m)  and  depends  upon  atmospheric  turbulence  and  stability,  and  surface  characteristics 
(Cheng  et  al.,  2001).  The  atmospheric  resistance  term,  Ra,  is  used  to  parameterize  the 
rate  of  pollutant  transfer  within  this  sub-layer  as  a  function  of  atmospheric  turbulence  and 
stability,  and  surface  characteristics  (Wesely  and  Hicks,  1977). 


Figure  1.       Relative  locations  where  dry  deposition  resistance  factors  Ra,  Rb?  and 
Rc  apply. 


Boundary-layer  Resistance  (Rb).  The  boundary  layer  is  a  thin,  non-turbulent  layer  that 
develops  just  above  the  surface.  The  depth  of  this  layer  is  in  terms  of  millimeters  (mm). 
For  rough  surfaces,  this  layer  is  constantly  changing  and  Hicks  (1982)  reported  that  is 
likely  to  be  intermittently  turbulent.  The  rate  of  pollutant  transfer  within  this  layer  is 
determined  by  molecular  diffusion  for  gases  and  Brownian  diffusion  and  inertial 
impaction  for  particles.  The  boundary-layer  resistance  term,  Rb,  is  usually  parameterized 
in  terms  of  the  Schmidt  number  (viscosity  of  air  divided  by  the  diffusivity  of  the 
pollutant)  and,  for  particles,  the  Stokes  number  (which  is  a  function  of  the  gravitation 
settling  velocity,  friction  velocity,  and  the  viscosity  of  air). 

Surface  Resistance  (Rc).  Vegetation  is  a  major  sink  for  soluble  or  reactive  gaseous 
pollutants.  After  passing  through  stomata  of  vegetation,  soluble  pollutants  dissolve  in 
moist  mesophyll  cells  in  the  interior  of  leaves  (Wesely  and  Hicks,  1977).  Reactive 
pollutants,  e.g.  ozone,  may  also  interact  with  the  exterior  (cuticle)  of  the  leaves.  Due  to 
response  of  the  stomata  to  external  factors  such  as  moisture  stress,  temperature,  and  solar 
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radiation,  resistance  in  the  vegetation  layer  can  exhibit  significant  diurnal  and  seasonal 
variability.  The  surface  resistance  term,  Rc,  is  usually  parameterized  in  terms  of  the  three 
main  pathways  for  uptake/reaction  of  a  pollutant  within  the  vegetation  or  surface  (Wesely 
and  Hicks,  1977): 

•  Transfer  through  stomatal  pores  and  dissolution  or  reaction  in  the  mesophyll  cells. 

•  Reaction  with  or  transfer  through  the  leaf  cuticle. 

•  Transfer  into  the  ground/water  surface. 

Cheng  et  al.  (2001)  provide  additional  information  on  components  that  are  measured 
and/or  observed  in  estimating  the  surface  resistance  factor. 

Atmospheric  Pollutants.  Atmospheric  species  commonly  measured  for  dry  deposition 
using  the  inference  method  are  presented  in  Table  1 .  Wesely  and  Hicks  (2000)  report 
that  NO  dry  deposition  is  usually  negligible  because  of  its  low  solubility  and  low 
oxidizing  capacity.  It  is  usually  not  considered  for  measurement.  Cheng  et  al.  (2001) 
recommend  that  gaseous  ammonia  (NH3)  not  be  considered  when  estimating  dry 
deposition  because  sufficient  understanding  of  its  biochemistry  has  yet  to  be  achieved. 


Table  1.        Species  commonly  measured  for  dry  deposition  using  the  inference 
 method.  

•  Sulphur  species  (gaseous  SO2  and  S04^"  in  particulate  matter). 

•  Nitrogen  species  (gaseous  NO2,  acidic  gases  HNO3  and  HNO2,  and  particulate  NH4'^  and 
NO3-). 

•  Particulate  base  cations  (Na^,  Mg^^,  Ca^^  and  K^).  Co-deposition  of  these  base  cations 
results  in  a  reduction  of  the  amount  of  deposited  acidity. 


Particle-associated  species  derived  from  erosion  of  soil  or  plant  material  (Na^,  K"^,  Mg^^, 
and  Ca^"^)  tend  to  reside  on  larger  airborne  particles  (e.g.  >2  pm)  (Lovett,  1994).  The 
majority  of  airborne  mass  of  NH/,  S04^',  and  H^  reside  on  submicrometer  aerosols. 
Thus  collecting  PMio  or  larger-sized  airborne  particles  (TSP)  provide  more  efficient 
capture  of  particle-associated  species  derived  from  erosion  of  soil  or  plant  material. 

Concentrations  of  the  eleven  species  are  combined  into  Equation  2  to  estimate  potential 
dry  acid  input  surface  load  (PAIdry)  in  kilogram  hydrogen  equivalents  (Cheng  et  al., 
2001): 

PAL  =     [^^2]  ,  M  ,  [hNO,]  ^  [HNO,]  ^  ^  [SOI]  ^  [nO^]  ^  [nh:] 
64         46  47  63  96         62  18 


39       23  40  24 


[kg  H^/ha/yr] 
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Meteorological  Parameters.  Continuous  measurement  of  numerous  meteorological 
variables  is  necessary  to  allow  estimation  of  dry  deposition  of  gaseous  and  particulate- 
bound  species  to  specific  surfaces  using  the  inference  method.  Meteorological  variables 
ultimately  required  are  one-hour  or  15-minute  standard  deviation  of  wind  direction,  wind 
speed,  solar  radiation,  and  air  temperature  at  standard  height  (10  m)  and  near  the  surface 
(2  m)  (EPCM,  2000).  These  temperatures  are  used  to  establish  atmospheric  stability. 
The  presence  or  absence  of  a  wet  surface  also  affects  dry  deposition.  Consequently, 
surface  wetness  sensing  and  relative  humidity  measurements  are  required.  In  summary, 
continuous  measurement  of  the  following  meteorological  variables  is  required  for 
estimating  gaseous  and  particulate-bound  species  dry  deposition  using  the  inference 
method: 

•  Wind  speed  and  wind  speed  standard  deviation 

•  Wind  direction  and  wind  direction  standard  deviation 

•  Solar  radiation 

•  Relative  humidity 

•  Surface  wetness 

•  Air  temperature  at  standard  height  (10  m) 

•  Difference  in  air  temperature  at  standard  height  and  surface  (taken  as  2  m  above 
ground). 

2.2      Methods  Description 

2.2.1  Alberta  En  vironment  Method 

Alberta  Environment  uses  an  inferential  method  to  estimate  aerodynamic  (Ra),  boundary- 
layer  (Rb),  and  surface  (canopy)  resistances  (Rc)  (Cheng  et  al.,  2001).  The  method 
parameterizes  dry  deposition  in  terms  of  a  deposition  velocity  (Vd)  that  is  inferred  from 
measured  meteorological  parameters  and  surface  characteristics  (Equation  1). 

Hourly  deposition  fluxes  for  each  species  are  calculated  as  a  product  of  the  hourly  Vd  and 
the  corresponding  hourly  concentration  (C).  Total  daily  dry  deposition  is  the  sum  of 
deposition  over  twenty-four  hours.  In  turn,  monthly,  seasonal,  and  annual  totals  are 
derived  by  summing  all  daily  totals  for  a  month,  season  or  year,  respectively.  WBK 
(2008)  used  this  method  to  estimate  annual  sulphur  and  nitrogen  species  deposition  using 
datasets  at  the  Alberta  Environment  acid  deposition  monitoring  site  at  Beaverlodge, 
Alberta  for  the  1998  to  2002  period.  Parameters  and  equations  used  in  the  Alberta 
Environment  method  are  presented  in  Appendix  I. 

2.2.2  En  vironment  Canada  Method 

Environment  Canada  uses  another  variation  of  the  inferential  method  to  estimate  dry 
deposition  velocities  of  acidic  species  at  sites  in  their  national  dry  deposition  monitoring 
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network  (Zhang  et  al.,  2003a,  2003b,  2002a,  2002b,  2001a,  and  2001b).  Their  approach 
is  known  as  A  Unified  Regional  Air  Quahty  Modehng  System  (AURAMS)  (Zhang  et  al., 
2002a).  The  network  is  referred  to  as  the  Canadian  Air  and  Precipitation  Monitoring 
Network  (CAPMoN)  (http://www.msc-smc.ec.gc.ca/capmon/index_e.cfm). 

A  Big  Leaf  model  is  used  for  gaseous  dry  deposition  in  AURAMS  (Zhang  et  al.,  2002a). 
A  size-segregated  particle  dry  deposition  module  originally  developed  by  the 
Meteorological  Service  of  Canada  is  used  for  particle  dry  deposition  in  AURAMS 
(Zhang  et  al.,  2001). 

AURAMS  was  recently  revised  to  include  non-stomatal  resistance  parameterizations 
(Zhang  et  al.,  2003b).  The  Big  Leaf  model  developed  by  Zhang  et  al.  (2002a)  was 
developed  for  calculating  dry  deposition  velocities  for  more  than  40  gaseous  species  for 
AURAMS,  but  it  only  included  seasonally- adjusted  values  for  non-stomatal  resistance. 
The  revised  model  incorporates  non-stomatal  resistance  parameterizations  (Zhang  et  al., 
2003a;  Zhang  et  al.,  2002b).  Other  improvements  included  more  practical  treatment  of 
cuticle  and  ground  resistance  in  winter  and  handling  of  seasonally-dependent  input 
parameters. 

Similar  to  the  Alberta  Environment  method,  AURAMS  uses  1-hour  average  values  of 
meteorological  observations  for  estimating  resistance  terms.  For  a  typical  31 -day  month, 
31  X  24  =  744  different  hourly  meteorological  observations  are  used  to  compute  a  similar 
number  of  hourly  average  deposition  velocities  and  deposition  loadings  for  each  acidic 
species.  A  monthly  deposition  load  is  computed  by  sunaming  individual  hourly  average 
loadings  (deposition  velocity  x  concentration  for  each  hour).  Parameters  and  equations 
used  in  the  Environment  Canada  approach  are  presented  in  Appendix  n. 

2,2.3    Wood  Buffalo  En  vironmental  Association  Metliod 

EPCM  (2002,  2000)  developed  a  method  for  estimating  SO2  and  NO2  dry  deposition  for 
the  Wood  Buffalo  Environmental  Association  Terrestrial  Environmental  Effects 
Monitoring  (TEEM)  Committee.  Aerodynamic  (Ra)  and  boundary  layer  (Rb)  resistances 
are  determined  using  simplified  empirical  relationships.  Surface  (canopy)  resistance  (Rc) 
is  determined  using  an  approach  very  similar  to  that  used  by  Alberta  Environment.  These 
relationships  are  described  in  Appendix  EL 

The  method  uses  fifteen-minute  average  values  of  meteorological  observations  to 
estimate  resistance  terms.  For  a  typical  31-day  month,  31  x  24  x  4  =  2,976  different  15- 
minute  meteorological  observations  are  used  to  compute  a  similar  number  of  15-minute 
average  deposition  velocities  and  deposition  loadings  for  each  acidic  parameter.  A 
monthly  deposition  load  is  computed  by  summing  individual  15-minute  average  loadings 
for  that  month. 
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2.3      Methods  Assumptions 

WBEA  data  gaseous  and  particulate  species  and  meteorological  data  at  the  Fort  McKay 
air  monitoring  station  are  shown  in  Table  2.  The  location  of  AMS  #1  relative  to  oil  sand 
lease  hold  areas  in  WBEA  zone  is  shown  in  Figure  2.  Hourly  data  for  parameters  listed 
in  Table  2  were  received  in  electronic  form  from  WBEA.  These  data  were  compiled, 
processed,  and  evaluated  using  the  three  inferential  methods  mentioned  previously. 

A  number  of  assumptions  were  required  beyond  what  each  method  offered  as  default 
assumptions  or  boundary  conditions  in  order  to  enable  modeling  to  be  performed.  These 
assumptions  are  listed  below.  The  intent  for  presenting  these  assumptions  is  to  allow 
others  to  conduct  similar  model  evaluations  using  the  same  (or  other)  inferential  methods 
and  dataset. 

1 .       Default  characteristics  (specific  location  and  land  use  category)  of  the  Fort 
McKay  air  monitoring  station  (AMS  #1): 

Latitude  57.2       South  Negative 

Longitude  -111.6     West  Negative 

Time  zone  -5  West  Negative 

Reference  height  (z)  10  m 

Land  Use  Category  (LUC)  Coniferous  Forest 


Table  2.        Dry  deposition  model  parameters  for  Wood  Buffalo  Environmental 
Association  Fort  McKay  air  monitoring  station  -  AMS  #1  -  for  the 
period  1  January  2003  to  31  December  2003. 


Components 

Parameters 

Combustion  gases 

Continuous  (hourly)  SO2  and  NO2 

Acid  gases 

24-hour  integrated  annular  denuder  sample  every  6'*^  day  for  HNO3,  HNO2,  and 

NH3 

Particulate-associated 

24-hour  integrated  PM2.5  and  PMio  sample  every  6""  day  for  NH^^.  SO4" ,  and 

acidic  species 

NO3 

Particulate-associated 

24-hour  integrated  PM25  and  PMio  sample  every  6"^  day  for  Na"^,  K"^,  Mg""^,  and 

base  cation  species 

Ca^" 

Meteorological 

Wind  speed  and  wind  speed  standard  deviation 

parameters 

Wind  direction  and  wind  direction  standard  deviation 

(recorded  hourly) 

Solar  radiation 

Relative  humidity 

Surface  wetness 

Air  temperature  at  standard  height  (10m) 

Difference  in  air  temperature  at  standard  height  and  surface  (taken  as  2  m  above 

ground) 
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2.        Data  preprocessing  procedures: 


a.  Invalid  or  missing  hourly  data  values  for  were  replaced  by  the  month's 
hourly  median  value  for  two  or  more  successive  invaUd  or  missing  hourly 
data  points. 

b.  Invalid  or  missing  hourly  data  values  were  replaced  by  an  average  of 
previous  and  succeeding  hourly  value  for  one  invalid  or  missing  hourly  data 
point. 

c.  Intermittent  hourly  data  (i.e.  samples  collected  for  24  hours  every  6**^  day) 
were  assumed  to  have  the  same  value  for  whole  6-day  sampling  period. 

d.  Hourly  NO2  and  SO2  concentrations  were  converted  from  ppbv  to  |Xg/m^  at 
standard  conditions  (25  °C  and  101.325  kPa). 

3.  Solar  and  earth  system  coordinates: 

Solar  zenith  angle  (0)  and  solar  radiation  (SR)  were  used  in  the  Environment 
Canada  (WBEA)  model  and  Wood  Buffalo  Environmental  Association  (WBEA) 
model.  A  Visual  Basic  program  was  developed  to  compute  those  parameters. 
The  results  of  this  program  were  confirmed  with  data  obtained  from  the 
Astronomical  Applications  Department,  US  Naval  Observatory 
(http://aa.usno.navy.mil/). 

In  general,  if  no  solar  radiation  data  were  available,  solar  radiation  (SR)  was 
assumed  equal  to  1000*cosine  of  solar  zenith  angle  (i.e.  1000*cos^. 

4.  Alberta  Environment  model  general  assumptions: 

Model's  equations  were  adapted  from  Cheng  et  al.  (2005)  and  are  listed  in 
Appendix  I.  These  equations  were  programmed  in  MS  EXCEL®  using  Visual 
Basic. 

a.  Surface  resistance  (Rc)  for  all  species  and  boundary-layer  resistance  (Rb)  for 
particles  were  weighted  based  on  54°N  latitude  for  each  season. 

b.  HNO3  was  assumed  to  have  the  same  Rb  as  SO2. 

c.  NH4^  was  assumed  to  have  the  same  Rb  as  SO/'. 

d.  Rb  for  NO3",  Na^,  K^,  Ca^^,  and  Mg^^  assumed  to  be  the  same  and  half  of  the 
Rb  of  SO4'". 

e.  Rc  for  SO2  and  NO2  were  based  on  default  values  listed  in  Table  4.4  and  4.5 
of  Cheng  et  al.  (2001). 

f.  Rc  for  HNO3  and  HNO2  was  assumed  to  be  equal  to  0. 1  s/cm  for  all  surfaces 
and  all  seasons. 

g.  Rc  for  all  particulate  species  assumed  equal  to  zero. 
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Figure  2.       Location  of  acid  deposition  monitoring  site  (Fort  McKay,  AMS  #1) 
among  oil  sand  lease  holders  in  Wood  Buffalo  Environmental 
Association  zone  (after  COST,  2006). 


h.  Surface  roughness  length  (zo)  was  computed  from  the  wind  speed  at 
reference  height  (u)  and  standard  deviation  of  the  wind  direction  (a^)  when 
wind  speed  >6  m/s. 

i.  No  available  conditions  were  available  for  zq  for  the  month  of  February, 
2003.  Therefore,  the  average  zq  of  January  and  March  was  used  to  represent 
conditions  for  this  month. 

j.     Boundary  conditions  (cut-off  points  used)  included: 

i.  If  horizontal  wind  speed  at  reference  height  (ii)  <0.36  km/h,  then  u  =  0 
and  deposition  velocity  (V^)  =  0. 

ii.  If  temperature  difference  (Td)  =  0,  then     =  0. 

iii.  Weighted-average  Rc  was  set  constant  for  each  season, 
k.     Surface  wetness  conditions: 

i.  If  precipitation  (Prec)  >  0. 1  mm/hour,  it  is  raining  and     =  (wet). 

ii.  If  relative  humidity  (RH)  >S1A%,  it  is  raining  and  Vd  =  (wet). 
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iii.    Vd  =  Wetness-weighted  V^. 

5.        Environment  Canada  model  general  assumptions: 

An  electronic  file  originally  programmed  using  FORTRAN  was  reprogranmied  in 
MS  EXCEL®  using  Visual  Basic.  Details  of  model  equations  are  listed  in 
Appendix  II. 

Gaseous  model  component  assumptions: 

a.  Daily  precipitation  amounts  for  the  Environment  Canada  Fort  McMurray 
Airport  meteorological  station  were  obtained  from  Environment  Canada 
(Edmonton,  AB).  The  squared  hourly  cloud  fraction  was  used  as  an  hourly 
weight  to  convert  daily  precipitation  to  hourly  precipitation. 

b.  Boundary  conditions  (cut-off  points  used)  included: 

i.  If  w  <3.6  km/h,  set  u  =  3.6  km/h. 

ii.  If  r^<lE"^°,  setr^=  lE-^°. 

iii.  If  Bulk  Richardson  number  (Ri)  >0  and  solar  radiation  (SR)  >0  w/m^, 
SQtRi=lE-^\ 

iv.  If  Monin-Obukhov  Length  (L)  >5,  set  L  =  5. 
V.     If  aerodynamic  resistance  (Ra)  <5,  set  Ra  =  5. 

vi.  If  Ra  >  1 000,  set  Ra  =  1 000. 

vii.  If  Rc<10,  setRc=  10. 

c.  Surface  wetness  conditions: 

i.  If  precipitation  (Free)  >0.1  mm/h,  it  is  raining. 

ii.  If  RH  >dew  point,  set  RH  to  dew  point. 

iii.  If  rainy  or  dew-point  condition  exist,  fraction  of  stomatal  blocking 
under  wet  conditions  (Wsd  =  (Srad  -  200)/800. 

iv.  If  Wst  <0.5,  set  Wst  =  0.5. 

V.     Wyrweighted  Rc  used  to  calculate  Vj. 

Particulate  model  component  assumptions: 

a.  Particulate  matter  size  was  assumed  to  follow  lognormal  distribution  from 
0.001  nmto  1  0|Lim. 

b.  S04^'  in  particulate  form  was  assumed  to  follow  lognormal  distribution  with 
MMD  =  0.35  |am  and  GSD  =  2.0. 

c.  NH4^  was  assumed  to  have  same  Vd  as  S04^". 

d.  Na^  in  particulate  form  was  assumed  to  follow  lognormal  distribution  with 
MMD  =  5.12  |nm  and  GSD  =  2.64. 

e.  NO3 ,  K^,  Ca^^,  and  Mg^^  were  assumed  to  have  same  Vd  as  Na"^. 

f.  Boundary  conditions  (cut-off  points  used)  included: 
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i.  If  u  <3.6  km/h,  set  U  =  3.6  km/h. 

ii.  If  r^<lE  '",  setr,/=  IE 

iii.  If  Ri  >0  and  SR  >0  w/m",  set  Ri  =  1 E  '  \ 

iv.  If  Monin-Obukhov  Length  (L)  >5,  set  L  =  5. 

V.     If  aerodynamic  resistance  (Ra)  <5,  set  Ra  =  5.  p 

vi.  If  Ra  >999.9,  set  Ra  =  999.9. 
g.     No  criteria  are  used  for  surface  wetness  conditions. 

6.        Wood  Buffalo  Environmental  Association  (WBEA)  model: 

Equations  for  this  method  were  obtained  from  WBEA  Excel  files  and  related 
documentation  (EPCM  2002,  2000)  received  from  the  Wood  Buffalo 
Environmental  Association.  These  equations  are  listed  in  Appendix  III.  These 
equations  were  programmed  in  MS  EXCEL  using  Visual  Basic. 

This  model  was  developed  to  process  input  parameters  as  15-minute  averages. 
However  model  runs  used  I -hour  average  values  for  input  parameters  as  15- 
minute  average  values  were  unavailable. 

a.  Atmospheric  conditions: 

a.  stable  for  night  time. 

b.  unstable  for  day  time. 

b.  HNO3  assumed  to  have  same  Rb  as  SO2. 

c.  NH4^  assumed  to  have  same  Rb  as  S04^". 

d.  Rb  for  NO3",  Na^,  K^,  Ca^^,  and  Mg""^  assumed  to  be  the  same  and  half  of  Rb 
ofSO/' 

e.  Rc  for  SO2  and  NO2  were  based  on  default  values  listed  in  Table  4.4  and  4.5 
of  Cheng  et  al.  (2001). 

f.  Rc  for  HNO3  and  HNO2  assumed  equal  to  0. 1  s/cm  for  all  surfaces  and  all 
seasons. 

g.  Rc  for  all  particles  assumed  equal  to  zero. 

h.  Boundary  conditions  (cut-off  points  used)  included: 

a.  If  u  <0.36  km/h:  set  m  =  0 

b.  Ra<1000 

c.  Rb  <1000 

i.  Surface  wetness  conditions: 

a.  If  precipitation  (Free)  >0.1  mm/hour,  it  is  raining  and  Vd  =  (wet). 

b.  If  RH  >87.4%,  it  is  raining  and  Vd  =  Vd  (wet). 
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2.4 


Results  and  Discussion 


2.4,1  ModelResuIts 

Individual  Species  Deposition  -  Estimated  monthly  dry  deposition  rates  for  a  number  of 
sulphur  and  nitrogen  species  at  the  Fort  McKay  site  are  presented  for  each  inference 
method  in  Figures  3  to  9.  Monthly  dry  deposition  rates  for  sulphur  species  -  SO2  and 
SO4""  -  are  presented  in  Figures  3  and  4,  respectively.  Monthly  dry  deposition  rates  for 
nitrogen  species  -  NO2,  HNO3,  HNO2,  NOs',  and  NH/  -  are  presented  in  Figures  5  to  9, 
respectively.  All  deposition  rates  are  expressed  in  "kg  species/ha/month." 

As  a  general  observation,  all  three  inference  methods  show  a  similar  ability  to  represent 
changes  (variation)  in  dry  deposition  rates  from  month  to  month.  For  example,  Figure  3 
illustrates  monthly  SO2  dry  deposition  rates  as  kg  S02/ha/month.  The  4*,  7^^,  and  8^^ 
months  of  the  year  clearly  show  higher  deposition  rates  for  each  method  relative  to  other 
months,  hi  this  particular  instance  the  Environment  Canada  (ENVC)  method  yields 
higher  monthly  SO2  deposition  rates  relative  to  the  other  two  methods. 

Another  general  observation  is  made  is  that  the  Wood  Buffalo  Environmental 
Association  (WBEA)  method  yields  higher  monthly  deposition  rates  for  most  species 
modeled  using  the  2003  dataset.  This  includes  S04^'  (Figure  4),  NO2  (Figure  5),  HNO3 
(Figure  6),  HNO2  (Figure  7),  and  NH/  (Figure  9). 

Potential  Acid  Input  -  Estimated  monthly  total  potential  acid  input  (PAI)  associated 
with  dry  species  at  the  Fort  McKay  site  are  presented  for  each  inference  method  in  Figure 
10.  These  rates  are  expressed  in  units  of  kg  H^/ha/yr  and  they  were  computed  taking  into 
account  all  sulphur,  nitrogen,  and  base  cation  species  as  shown  in  Equation  2.  Similar  to 
findings  reported  above.  Figure  10  shows  that  the  WBEA  method  gives  higher  PAI  rates 
associated  with  dry  species  during  the  winter  and  spring  months  using  the  2003  dataset. 

The  WBEA  method  gives  similar  rates  as  the  ENVC  and  AENV  methods  during  other 
seasons  (i.e.  summer  and  fall).  It  is  initially  suspected  that  a  simpler  computational 
approach  for  the  aerodynamic  (Ra)  and  surface  (Rb)  resistance  terms  in  the  WBEA 
method  (Appendix  III)  partially  accounts  for  this  departure  during  winter  and  spring 
months.  However,  further  testing  would  be  necessary  to  better  understand  whether  this  is 
the  case. 

A  general  observation  is  made  regarding  the  performance  of  each  method  and  their 
ability  to  represent  change  (variation)  in  monthly  total  PAI  rates  associated  with  dry 
species  shown  in  Figure  10.  The  ENVC  and  AENV  method  yield  comparable  monthly 
total  PAI  rates  associated  with  dry  species,  except  for  the  1 1^^  month  (November)  where 
the  AENV  yields  a  net  negative  rate.  The  exact  reasons  for  this  departure  in  the  1 1^^ 
month  are  unknown,  however  it  may  be  due  to  greater  equivalent  total  base  cation  species 
deposition  flux  relative  to  equivalent  total  sulphur  and  nitrogen  species  flux  during  the 
month. 
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Figure  3.  Monthly  sulphur  dioxide  (SO2)  dry  deposition  for  the  period  January 
to  December  2003  at  WBEA  Fort  McKay  air  monitoring  station  AMS 
#1  (expressed  in  units  of  kg  SOi/ha/month). 


(ENVC  =  Environment  Canada  method;  AENV  =  Alberta  Environment  method; 
WBEA  =  Wood  Buffalo  Environmental  Association  method) 
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Figure  4.       Monthly  sulphate  (SO4  ")  dry  deposition  for  the  period  January  to 

December  2003  at  WBEA  Fort  McKay  air  monitoring  station  AMS  #1 
(expressed  in  units  of  kg  S04^7ha/month). 


(ENVC  =  Environment  Canada  method;  AENV  =  Alberta  Environment  method; 
WBEA  =  Wood  Buffalo  Environmental  Association  method) 
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Figure  5.  Monthly  nitrogen  dioxide  (NO2)  dry  deposition  for  the  period  January 
to  December  2003  at  WBEA  Fort  McKay  air  monitoring  station  AMS 
#1  (expressed  in  units  of  kg  NOi/ha/month). 

(ENVC  =  Environment  Canada  method;  AENV  =  Alberta  Environment  method; 
WBEA  =  Wood  Buffalo  Environmental  Association  method) 
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Figure  6.       Monthly  nitric  acid  (HNO3)  dry  deposition  for  the  period  January  to 
December  2003  at  WBEA  Fort  McKay  air  monitoring  station  AMS  #1 
(expressed  in  units  of  kg  HNOa/ha/month). 

(ENVC  =  Environment  Canada  method;  AENV  =  Alberta  Environment  method; 
WBEA  =  Wood  Buffalo  Environmental  Association  method) 
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Figure  7.  Monthly  nitrous  acid  (HNO2)  dry  deposition  for  the  period  January  to 
December  2003  at  WBEA  Fort  McKay  air  monitoring  station  AMS  #1 
(expressed  in  units  of  kg  HNOi/ha/month). 

(ENVC  =  Environment  Canada  method;  AENV  =  Alberta  Environment  method; 
WBEA  =  Wood  Buffalo  Environmental  Association  method) 
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Figure  8.       Monthly  nitrate  (NO3")  dry  deposition  for  the  period  January  to 

December  2003  at  WBEA  Fort  McKay  air  monitoring  station  AMS  #1 
(expressed  in  units  of  kg  NOa'/ha/month). 


(ENVC  =  Environment  Canada  method;  AENV  =  Alberta  Environment  method; 
WBEA  =  Wood  Buffalo  Environmental  Association  method) 
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Figure  9.       Monthly  ammonium  (NH4^)  dry  deposition  for  the  period  January  to 


December  2003  at  WBEA  Fort  McKay  air  monitoring  station  AMS  #1 
(expressed  in  units  of  kg  NH4^/ha/month). 

(ENVC  =  Environment  Canada  method;  AENV  =  Alberta  Environment  method; 
WBEA  =  Wood  Buffalo  Environmental  Association  method) 
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Figure  10.      Monthly  potential  acid  input  (PAI)  values  associated  with  dry  species 
for  the  period  January  to  December  2003  at  WBEA  Fort  McKay  air 
monitoring  station  AMS  #1  (expressed  in  units  of  kg  H^/ha/month). 

(ENVC  =  Environment  Canada  method;  AENV  =  Alberta  Environment  method; 
WBEA  =  Wood  Buffalo  Environmental  Association  method) 
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It  is  observed  that  ENVC  model  parameters  have  more  boundary  conditions  (i.e.  cut-off 
points  and  lower/upper  hmits)  associated  with  their  use  (Appendix  II)  compared  to  the 
AENV  model  (Appendix  I).  Although  not  shown,  it  was  observed  that  the  AENV  model 
had  much  greater  variabihty  associated  with  hourly  model  outputs.  This  may  partially 
explain  a  net  negative  total  PAI  rate  during  the  1 1'*^  month  for  the  2003  dataset.  At  this 
time  the  influence  of  assigning  additional  boundary  conditions  associated  with  AENV 
model  input  parameters  on  variability  associated  with  model  outputs  is  unknown.  Further 
testing  would  be  required  to  better  understand  whether  this  has  an  ability  to  improve 
AENV  model  performance  relative  to  the  ENVC  model. 

Annual  total  PAI  rates  associated  with  dry  species  at  the  Fort  McKay  site  using  the  2003 
dataset  were  0.15  kg  HVha/yr  (ENVC  method),  0.14  kg  H^/ha/yr  (AENV  method),  and 
0.20  kg  H^/ha/yr  (WBEA  method).  The  higher  annual  rate  associated  with  dry  species 
observed  for  the  WBEA  method  relative  to  the  other  methods  is  entirely  consistent  with 
higher  individual  species  deposition  rates  observed  in  Figures  4  to  9. 

Limited  data  exist  to  enable  comparisons  of  these  results  with  other  observations. 
Environment  Canada  (2004)  reported  on  modeling  activities  in  the  2004  Alberta-wide 
acid  deposition  assessment.  Annual  average  concentrations  and  annual  depositions  of 
sulphur  species  were  calculated  for  a  1  °  by  1  °  latitude-longitude  grid  area  where  the  Fort 
McKay  site  is  located  using  the  REgional  Lagrangian  Acid  Deposition  (RELAD)  model. 
Modeling  was  undertaken  for  the  1971-2000  period  using  30  years  of  meteorological  data 
and  1995  emissions. 

Environment  Canada  (2004)  presented  data  that  indicate  a  1971  to  2000  (30-yr)  mean 
annual  sulphur  species  dry  deposition  rate  range  of  0.06  to  0.08  kg  H^/ha/yr  for  a  1°  by  1° 
latitude-longitude  grid  area  where  the  Fort  McKay  site  is  located.  Although  these  results 
are  not  directly  comparable,  annual  total  sulphur  species  dry  deposition  rates  at  the  Fort 
McKay  site  using  the  2003  dataset  were  0.053  kg  H^/ha/yr  (ENVC  method),  0.044  kg 
H'-Zha/yr  (AENV  method),  and  0.060  kg  H^/ha/yr  (WBEA  method). 

RWDI  (2004)  reported  historical  measured  and  predicted  NO2  dry  deposition  rates  for  the 
Fort  McKay  site.  RWDI  (2004)  reported  an  average  measured  NO2  dry  deposition  rate  of 
0.17  kg  H'-Zha/yr  (2.33  kg  N/ha/yr)  for  the  1999  to  2000  period.  RWDI  (2004)  also 
reported  an  average  modeled  NO2  dry  deposition  rate  of  0.14  kg  H^/ha/yr  (1.93  kg 
N/ha/yr)  for  the  1998  to  2002  period.  Annual  total  NO2  dry  deposition  rates  at  the  Fort 
McKay  site  using  the  2003  dataset  were  0.10  kg  H^/ha/yr  (ENVC  method),  0.08  kg 
H^/ha/yr  (AENV  method),  and  0. 1 2  kg  H^/ha/yr  (WBEA  method). 

2,4,2    Comparison  of  Inference  Methods 

A  general  evaluation  of  model  complexity  and  performance  was  recorded  after  work  was 
completed  with  each  inferential  method  for  estimating  dry  deposition.  A  relative 
comparison  of  selective  features  of  the  three  inference  methods  is  presented  in  Table  3. 
Overall,  it  is  the  view  of  the  report  authors  that  the  Environment  Canada  (ENVC)  method 
is  more  sophisticated  (complicated)  given  that  it  has  the  greatest  number  of  input 
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parameter  requirements.  Although  not  shown,  it  also  required  making  a  number  of 
additional  assumptions  for  some  of  the  input  parameters  before  modeling  could  be 
performed. 


Table  3.        Relative  comparison  of  selective  features  of  three  inference  methods 
for  estimating  dry  deposition. 


ENVC 

AENV 

WBEA 

Number  of  Parameters 

18 

9 

7 

Relative  level  of  complexity 

High 

Medium 

Low 

Relative  computing  time 

High 

Low 

Low 

requirements 

Method  for  calculating  Ra 

Resistance  model 

Resistance  model 

Empirical  equation 

Method  for  calculating  Rb 

Resistance  model 

Resistance  model  for 

Resistance  model  for 

gas  and  empirical 

gas  and  empirical 

constant  for  particles 

constant  for  particles 

Method  for  calculating  Rc 

Resistance  model 

Empirical  constant 

Empirical  constant 

Major  assumptions  required 

Meteorological  data: 

n/a 

n/a 

houriy  precipitation 

obtained  from  daily 

totals  in  vague  way 

Particle  mass 

distribution  was 

assumed 

Number  of  land  use  categories 

26 

8 

8 

(LUCs) 

Seasons 

Function  of  latitude 

Function  of  month, 

Function  of  month 

and  month 

4  seasons  total 

Number  of  gas  species 

31 

4 

4 

Assumptions  for  particulate 

Based  on  distribution 

7 

7 

species 

assumptions 

Wetness  correction 

Corrected  based  on 

Time  weighted  by 

Corrected  by 

rain  and  dew 

combined  surface 

combined  surface 

condition,  respectively 

wetness  and  related 

wetness  and  related 

humidity  conditions 

humidity  conditions 

Basis  for  leaf  area  index  (LAI) 

Function  of  LUC  and 

n/a 

n/a 

season 

Basis  for  surface  roughness 

Function  of  LUC  and 

Function  of  wind 

n/a 

Season 

speed  and  wind 

direction  deviation 

Basis  for  stability  condition 

Temperature  difference 

Temperature 

Day  vs.  night 

with  temperature 

difference 

(unstable  vs.  stable 

gradient  correction 

assumption) 

ENVC  =  Environment  Canada  method,  AENV  =  Alberta  Environment  method, 
WBEA  =  Wood  Buffalo  Environmental  Association  method 
n/a  =  not  applicable 
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The  Alberta  Environment  (AENV)  method  provided  comparative  results  to  the  ENVC 
method  based  on  the  2003  Fort  McKay  dataset.  This  method  is  less  complicated 
compared  to  the  ENVC  method,  and  it  requires  less  input  data  and  computing  time. 
Finally,  the  Wood  Buffalo  Environmental  Association  (WBEA)  method  is  simple  and 
straight-forward  to  program.  However  as  discussed  previously,  this  method  yielded 
higher  monthly  deposition  rates  for  most  species  modeled  using  the  2003  dataset. 

Although  a  Hmited  one-year  dataset  was  available  to  evaluate  these  models,  it  is  the  view 
of  the  report  authors  that  the  AENV  method  offers  a  practical  choice  for  routine 
computation  of  sulphur  and  nitrogen  species  dry  deposition  in  Alberta.  Further  testing 
using  different  datasets  elsewhere  in  Alberta  is  recommended  to  support  this.  If  the 
ENVC  method  can  be  shown  to  provide  comparative  results  to  the  AENV  method  using 
different  datasets  elsewhere  in  Alberta,  it  offers  a  logical  choice  for  periodic  use  as  a 
check  on  the  AENV  method. 

2,4.3  Relationships  of  Dry  Deposition  for  Sulpliur  and  Nitrogen  Species 

Contributions  of  SO2  and  NO2  deposition  in  total  sulphur  and  nitrogen  species  deposition 
at  the  Fort  McKay  site  were  examined  for  the  one-year  dataset.  Annual  sulphur  and 
nitrogen  species  deposition  data  were  analyzed  to  determine  the  ratio  of  SO2  deposition 
to  total  sulphur  species  deposition  and  of  NO2  deposition  to  total  nitrogen  species 
deposition.  These  results  are  presented  in  Table  4  (SO2  to  total  S  species  deposition 
ratio)  and  Table  5  (NO2  to  total  N  species  deposition  ratio),  respectively. 

Table  4.        Ratio  of  annual  SO2  deposition  to  total  sulphur  species  deposition  at 
WBEA  Fort  McKay  air  monitoring  station  during  2003. 


Method 

SO2 

S04^" 

S02/S,ot 

kg  S/ha/yr 

kg  S/ha/yr 

2.4.3.1  Ratio 

ENVC 

1.489 

0.101 

0.94 

AENV 

0.975 

0.222 

0.82 

WBEA 

1.197 

0.361 

0.77 

SO, 

Note: 

  Ratio  = 

so^  +  sol' 

ENVC 

=  Environment  Canada  method 

AENV 

=  Alberta  Environment  method 

WBEA 

=  Wood  Buffalo  Environmental  Association  method 

Results  presented  in  Table  4  and  5  using  the  Alberta  Environment  method  (AENV)  are 
discussed  here  because  previous  work  by  WBK  (2008)  used  the  same  method  to  estimate 
ratios  at  another  acid  deposition  monitoring  site  in  Alberta  -  Beaverlodge.  Estimates  of 
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annual  S  species  deposition  for  the  Fort  McKay  site  indicate  that  about  82%  is  in  the 
form  of  gaseous  SO2  with  the  remainder  as  particulate  sulphate  (18%). 


WBK  (2008)  reported  on  annual  deposition  loadings  of  S  and  N  species  during  1998  to 
2002  at  a  former  Alberta  Environment  Beaverlodge  acid  deposition  monitoring  site.  The 
station  was  located  west  of  Grande  Prairie  at  an  Agriculture  and  Agri-food  Canada 
Research  Farm.  Estimates  of  annual  S  species  deposition  (Table  6)  consistently  showed 
that  about  80%  was  in  the  form  of  gaseous  SO2  with  the  remainder  as  particulate 
sulphate.  These  data  compare  well  with  that  observed  at  Fort  McKay  during  2003  (82%). 

Table  5.        Ratio  of  annual  NO2  deposition  to  total  nitrogen  species  deposition  at 
WBEA  Fort  McKay  air  monitoring  station  during  2003. 


Method 

NO2 
Kg  N/ha/yr 

HNO3 
kg  N/ha/yr 

HNO2 
kg  N/ha/yr 

NOs' 
kg  N/ha/yr 

NH/ 
kg  N/ha/yr 

N02/Nto, 
2.4.3.2  Ratio 

ENVC 

1.368 

0.430 

0.083 

0.038 

0.112 

0.67 

AENV 

1.132 

0.524 

0.250 

0.051 

0.250 

0.51 

WBEA 

1.632 

0.850 

0.463 

0.082 

0.404 

0.48 

NO,  NO, 
Note:   Ratio  =  

N,^,  NO2  +  HNO,  +  HNO2  +         +  ^^^3' 

ENVC  =  Environment  Canada  method 
AENV  =  Alberta  Environment  method 
WBEA  =  Wood  Buffalo  Environmental  Association  method 


Estimates  of  annual  N  species  deposition  at  the  Fort  McKay  site  showed  that  about  5 1  % 
was  from  NO2  with  the  remainder  as  HNO3  and  HNO2  (35%)  and  particulate  ammonium 
and  nitrate  (<14%).  Estimates  of  annual  N  species  deposition  at  the  Beaverlodge  site 
(Table  7)  showed  that  about  35  to  50%  was  from  NOx  (NO  +  NO2)  with  the  remainder  as 
HNO3  and  HNO2  (-40  to  60%)  and  particulate  ammonium  and  nitrate  (<10%)  (WBK, 
2008).  Again,  results  for  the  Beaverlodge  site  on  the  same  order  as  that  observed  at  the 
Fort  McKay  site. 

WBK  (2008)  reported  on  work  of  Peake  and  Davidson  (1990)  related  to  estimated  annual 
dry  deposition  of  N  species  (NOx,  HNO2,  HNO3,  and  NO3")  in  the  south  western  region  of 
Alberta.  This  region  stretches  east  from  the  Great  Divide  of  the  Rocky  Mountains  to  the 
plains  of  southern  Alberta,  80  km  east  of  Calgary,  as  discussed  by  Peake  and  Davidson 
(1990).  These  estimates  were  based  upon  measurements  made  at  Crossfield  east  and 
west,  and  Fortress  Mountain  monitoring  sites  during  1985  to  1987  as  part  of  the  Alberta 
Government/Industry  Acid  Deposition  Research  Program  (ADRP).  Using  data  reported 
by  Peake  and  Davidson  (1990).  WBK  (2008)  estimated  that  about  32%  of  N  deposition 
was  from  NOx  (NO  +  NO2)  with  the  remainder  as  nitric  and  nitrous  acid  (-63%)  and 
particulate  nitrate  (-5%). 


Dry  Deposition  Monitoring  Method  in  Alberta 


25 


Table  6. 


Ratio  of  annual  SO2  deposition  to  total  sulphur  species  deposition  at 
Beaverlodge,  Alberta  (after  WBK,  2008). 


Year 

kg  S/ha/yr 

<50  ^' 

kg  S/ha/yr 

OU2/Otol 

2.4.3.3  Ratio 

1998 

0.368 

0.083 

0.82 

1999 

0.305 

0.072 

0.81 

2000 

0.294 

0.064 

0.82 

2001 

0.317 

0.069 

0.82 

2002 

0.317 

0.072 

0.81 

Note: 

  Ratio  = 

SO, 

so,  +  " 

Bytnerowicz  et  al.  (1999)  as  cited  in  Bytnerowicz  et  al.  (2005)  reported  that  HNO3 
typically  provides  more  than  60%  of  all  dry-deposited  N  species  in  mixed  conifer  forests 
of  the  Los  Angeles  Basin  mountain  range  of  California. 


Table  7.        Ratio  of  annual  NOx  deposition  to  total  nitrogen  species  deposition  at 
Beaverlodge,  Alberta  (after  WBK,  2008). 


Year 

NOx 
Kg  N/ha/yr 

HNO3 
kg  N/ha/yr 

HNO2 
kg  N/ha/yr 

NOs" 
kg  N/ha/yr 

NH4* 
kg  N/ha/yr 

NOx/Ntot 
2.4.3.4  Ratio 

1998 

0.617 

0.370 

0.097 

0.030 

0.077 

0.52 

1999 

0.556 

0.482 

0.051 

0.034 

0.060 

0.47 

2000 

0.551 

0.412 

0.064 

0.029 

0.054 

0.50 

2001 

0.572 

0.600 

0.098 

0.035 

0.050 

0.42 

2002 

0.679 

0.952 

0.229 

0.041 

0.040 

0.35 

NO. 

Note:  NO,/Ntot  Ratio  = 


NO^  +  HNO,  +  HNO,  +  NH;  +  NO^ 
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3.0    FINDINGS  AND  RECOMMENDATIONS 


1 .  Complete  sulphur  and  nitrogen  gaseous  and  particulate  species  and 
meteorological  data  for  a  one-year  period  (2003)  from  the  Wood  Buffalo 
Environmental  Association  Fort  McKay  air  monitoring  station  were  evaluated 
using  the  three  methods.  All  three  methods  showed  a  similar  ability  to  represent 
variation  in  dry  deposition  rates  for  individual  species  from  month  to  month. 

2.  The  WBEA  method  yielded  higher  monthly  deposition  rates  for  most  species 
modeled  (SO4" ,  NO2,  HNO3,  HNO2,  and  NH/)  compared  to  the  AENV  and 
ENVC  methods.  The  WBEA  method  also  yielded  higher  potential  acid  input 
(PAI)  values  associated  with  dry  species  compared  to  the  AENV  and  ENVC 
methods  during  winter  and  spring.  The  WBEA  method  gave  similar  PAI  values 
associated  with  dry  species  during  summer  and  fall.  A  simpler  computational 
approach  for  aerodynamic  and  surface  resistance  terms  in  the  WBEA  method  may 
partially  explain  differences  observed  during  winter  and  spring.  Further  testing 
would  be  required  to  better  understand  whether  this  is  the  case. 

3.  Annual  total  PAI  associated  with  dry  species  for  2003  was  0. 14  kg  H'^/ha/yr 
(AENV  method),  0.15  kg  R-'/ha/yr  (ENVC  method),  and  0.20  kg  H^/ha/yr 
(WBEA  method).  The  higher  annual  rate  associated  with  dry  species  observed 
for  the  WBEA  method  relative  to  other  methods  is  consistent  with  higher  monthly 
deposition  rates  for  individual  species.  While  the  AENV  and  ENVC  methods 
yielded  comparable  monthly  PAI  values  associated  with  dry  species,  it  is 
observed  that  ENVC  input  parameters  have  more  boundary  conditions  (i.e.  cut-off 
points  and  lower/upper  limits)  associated  with  their  use  compared  to  AENV  input 
parameters. 

4.  The  ENVC  method  is  more  complicated  given  that  it  has  the  greatest  number  of 
input  parameter  requirements.  A  number  of  additional  assumptions  for  some  of 
the  input  parameters  are  also  required  before  modeling  can  be  performed.  The 
AENV  method  provided  comparative  results  to  the  ENVC  method.  The  AENV 
method  is  less  complicated  compared  to  the  ENVC  method,  and  it  requires  less 
input  data  and  computing  time.  The  WBEA  method  is  simple  and  straight- 
forward to  use.  However,  this  method  yielded  higher  monthly  deposition  rates  for 
most  species  modeled  compared  to  the  other  methods  using  the  2003  dataset. 

5.  The  dataset  used  to  evaluate  the  methods  was  short  -  one  year  -  and  month  to 
month  variation  or  lack  of  variation  in  species  deposition  rates  observed  using 
each  method  does  not  necessarily  provide  an  indication  of  what  may  happen  in 
other  years.  Although  only  a  one-year  dataset  was  available  for  evaluation, 
results  support  that  the  AENV  method  offers  a  practical  choice  for  routine 
computation  of  sulphur  and  nitrogen  species  dry  deposition  in  Alberta.  Further 
testing  using  different  datasets  elsewhere  in  Alberta  is  recommended  to  support 
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this.  If  the  ENVC  method  can  be  shown  to  provide  comparable  resuhs  to  the 
AENV  method  using  different  datasets,  it  would  offer  a  logical  choice  for 
periodic  use  as  a  check  on  the  AENV  method. 

6.  Evaluation  of  contributions  of  SO2  and  NO2  deposition  in  total  sulphur  and 
nitrogen  species  deposition  at  the  Fort  McKay  site  indicated  that  about  827r  is  in 
the  form  of  gaseous  SO2  with  the  remainder  as  particulate  SO4"  (18%).  These 
results  are  comparable  to  annual  deposition  loadings  of  sulphur  species  during 
1998  to  2002  at  the  former  Alberta  Environment  Beaverlodge  acid  deposition 
monitoring  site  (80%  in  the  form  of  gaseous  SO2  and  20%  as  particulate  SO4" ). 

Estimates  of  annual  nitrogen  species  deposition  at  the  Fort  McKay  site  indicated 
that  -51%  was  from  NO2  with  the  remainder  as  HNO3  and  HNO2  (35%)  and 
particulate  NH/  and  NO3"  (<14%).  Estimates  of  annual  nitrogen  species 
deposition  at  the  Beaverlodge  site  during  1998  to  2002  showed  that  about  35  to 
50%  was  from  NO  and  NO2  with  the  remainder  as  HNO3  and  HNO2  (40  to  60%) 
and  particulate  NH/  and  NO3  (<10%). 

7.  Future  work  evaluating  these  inference  methods  should  consider  the  influence  of 
assigning  additional  boundary  conditions  to  AENV  model  input  parameters  on 
variability  associated  with  model  outputs.  Specifically,  assigning  boundary 
conditions  to  AENV  input  parameters  consistent  with  default  ENVC  model 
assumptions  should  be  investigated  to  better  understand  whether  these  changes 
have  the  ability  to  improve  AENV  model  performance  relative  to  the  ENVC 
model. 
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APPENDIX  I 

Alberta  Environment  Calculation  Methods  for  Gases  and 

Particulates 

(after  Cheng  et  al.,  2001) 


Parameters: 

PAI 

Potential  Acid  Input  (kg  H^/ha/yr) 

[X] 

concentration  of  X  chemical  species  deposited  (kg/ha/yr) 

F 

dry  deposition  flux  (ug/m^/s) 

Vd 

deposition  velocity  (m/s) 

c 

concentration  (ug/m^) 

Ra 

aerodynamic  resistance  (s/m) 

Rb 

boundary-layer  resistance  (s/m) 

Rc 

surface  resistance  (s/m) 

k 

von  Karman  constant  (0.4) 

* 

u 

friction  velocity  (m/s) 

z 

reference  height  (10  m) 

zo 

surface  roughness  length  (m) 

¥ 

integrated  stability  correction  term 

L 

Monin-Obukhov  length  scale 

u 

wind  speed  (m/s) 

Standard  deviation  of  wind  direction  (rad) 

Ri 

Bulk  Richardson  number 

a 
& 

pravitational  arrplpration  (^Q  Rl  m/s^^ 

Td 

temperature  difference  between  10  and  2  m  (Tio  -  T2) 

T2 

temperature  at  2  m  (Kelvin) 

H 

sensible  heat  flux  (w/m^) 

B 

defined  equation  described  herein 

dynamic  viscosity  of  air  (18.0  x  10"^  N-s/m*^  at  1  atm  and  25  °C) 

p 

density  of  air  (1.18  kg/m^  at  1  atm  and  25°C) 

D 

diffusion  coefficient  of  the  substance  of  interest  (cm^/s) 

Pr 

Prandtl  number  for  air  (0.72) 

Tl/pD 

Schmidt  number 

RH 

relative  humidity 

SW 

soil  wetness 
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Equations: 


_  [so,]  ^  [no,]  ^  [hNO,]  ^  [hNO,]  ^  ,[.SO;-]  ^  [nQ-]  ^  [nh;]  deposit 


ion 


PAIdry 

^       64  46  47  63  96  62  18  concentrations 

^  39        23  40  24 

F  =  VhC 


Vd  = 


Summary  of  Species  Specific  Deposition  Velocity  Formulae: 

1 


{SO 2 


^  (^^2 )  1^  Rc  is  treated  as  being  negligible  for  nitric  and 

nitrous  acid. 


diHNO^] 


V 


d[SOf,NH;  J?    _L  J? 


d[MO^Xa^-,Mr^,K^,Ma^}  +  (0.5  X .  , 

Aerodynamic  Resistance: 

Ra  =  —tI  ^^~~yn  ~  \  }      (Ra  is  infinite  and     =  0  when  u  and  Td  are  zero) 


kU     \      Zn  \L 


u  =         (this  relationship  is  used  as  an  initial  estimate  of  u  to  calculate  zq,  a  more 
precise  value  of  u  is  calculated  after  zq  is  obtained  -  refer  to  below) 
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OAu' 


zo  =  z  e  ^  (calculated  as  a  monthly  average  using  data  where  the  wind  speed  is  >6 
m/s) 


A  more  precise  value  of  u"^  calculated  after  zo  is  obtained  (based  on  atmospheric 
conditions): 


Stable  conditions       Rj  >  0       u*  =ku<\n 


r  \ 
z 


(1+4.7/?/; 


Unstable  conditions    Ri  <  0       w  * 


ku 


In 


r  \ 
z 


9ARi 


1  -t-  (7.45) 


Neutral  conditions     Rj  =  0  u* 


ku 


In 


z 


Calculation  of  \|/  (based  on  atmospheric  conditions): 


5z 


T^u 


Stable  conditions       vj/  =  -  —        where  L  -  — — 

L  kHg 


Unstable  conditions    v|/  =  2  ln< 


1  +  Jl- 


15z 


Neutral  conditions     \|/  =  0 
Calculation  of  H  (based  on  atmospheric  conditions): 


Neutral  and  Stable  conditions  H 


uT^ 
0.74 


In 


z 


(1  +  4.7/?/)' 
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2 


Unstable  conditions 


where  B  =  9.4 


In  ^/ 


0.74 


In 


9.4/?/ 

[1+5.35; 


Boundary  Layer  Resistance: 
Rb  (gases): 


Rb  -  Rh  - 


2 


77  1 

—  X  

^pD  Pr 


Rb  =  ^  for  SO2  and  HNO3 
u 


Rb=  forNOs 


Rb=  ^  forHN02 


Rb  (particulates): 

Rb  values  for  particulate  sulphate  are  obtained  from  scientific  literature  for  daytime  and 
nighttime  as  a  function  of  surface  type  and  weighted  according  to  average  day  length  for 
each  month  at  a  mid- Alberta  latitude  location  (54°N  latitude)  after  Cheng  and  Angle 
(1993)  as  cited  in  Cheng  et  al.  (2001). 

Boundary-Layer  Resistance  (s/cm)  for  Particulate  Sulphate,  Day  Length  Weighted 
Averages  at  54°N  Latitude  for  the  Middle  of  Each  Month.  


Winter 


Spring 


Summer 


Autumn 


Surface  Type 

(Dec,  Jan,  Feb) 

(Mar,  Apr,  May) 

(Jun,  July,  Aug) 

(Sep.  Oct,  Nov) 

Dry 

Wet 

Dry 

Wet 

Dry 

Wet 

Dry 

Wet 

Deciduous  Forest 

16.9 

0 

5.4 

0 

1.3 

0 

3.2 

0 

Coniferous  Forest 

2.5 

0 

2.7 

0 

1.9 

0 

2.3 

0 

Wetland/Swamp* 

20.4 

0 

3.8 

0 

2.6 

0 

3.2 

0 

Grassland* 

20.4 

0 

5.6 

0 

3.9 

0 

4.7 

0 

Cropland* 

20.4 

0 

9.0' 

0 

3.9 

0 

7.9* 

0 

Urban^ 

33.9 

0 

10.9 

0 

2.6 

0 

6.3 

0 
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Open  Water  0  0  0  0  0  0  0 

Snow/Ice  204  0  -  -  -  

*  in  winter,  wetland,  grassland,  and  cropland  treated  as  a  snow  surface.  '  bare  soil  and  active 

growth. 

^  bare  soil  and  senescent  growth.  ^  consists  of  a  mixture  of  deciduous  forest  and 

buildings. 


Day  length  Weighted  Seasonal  Average R,=  i?!^^k,J+  X-^-!!L^{r^„J 

Day  length  =  0.133{:os"'  (-  tan(55°  )x  idiXv{Solar  Declination)} 

360x(2M  + Julian  DayY 


day  length 


Solar  Declination  =  23.45<!  sin 


365 


Surface  Resistance: 


Bulk  surface  resistance  values  are  used  from  literature  as  a  function  of  surface  type, 
surface  wetness,  and  incident  radiation.  Day  length  weighted  average  Rc  values  for  SO2 
and  NO2  are  used  from  Voldner  et  al  (1986),  Arrit  et  al  (1987)  and  Walcek  et  al  (1986)  as 
cited  in  Cheng  et  al.  (2001): 

Day  Length  Weighted  Averages  Bulk  Surface  Resistance  (s/cm)  for  Sulphur  Dioxide 


(S02): 

Winter 

Spring 

Summer 

Autumn 

Surface  Type 

(Dec,  Jan,  Feb) 

(Mar,  Apr,  May) 

(Jun,  July,  Aug) 

(Sep,  Oct,  Nov) 

Dry 

Wet 

Dry 

Wet 

Dry 

Wet 

Dry 

Wet 

Deciduous  Forest 

10 

10 

4.7 

0 

3.5 

0 

7.9 

0.4 

Coniferous  Forest 

5 

5 

4.1 

0 

3.5 

0 

4.9 

0.2 

Wetland/Swamp* 

7 

1 

0.5 

0 

0.7 

0 

1 

0.1 

Grassland* 

7 

1 

1 

0 

1.3 

0 

2 

0.1 

Cropland* 

7 

1 

ot 

0 

2 

0 

2+ 

0.1 

Urban^ 

10 

2 

10 

0 

10 

0 

10 

0.1 

Open  Water 

0 

0 

0 

0 

0 

0 

0 

0 

Snow/Ice 

7 

1 

*  in  winter,  wetland,  grassland,  and  cropland  treated  as  a  snow  surface.  '  bare  soil  and  active 

growth. 

^  bare  soil  and  senescent  growth.  ^  consists  of  a  mixture  of  deciduous  forest  and 

buildings. 


Day  Length  Weighted  Averages  Bulk  Surface  Resistance  (s/cm)  for  Nitrogen  Dioxide 

(NO2):  

Winter  Spring  Summer  Autumn 

Surface  Type  (Dec,  Jan,  Feb)         (Mar,  Apr,  May)         (Jun,  July,  Aug)         (Sep,  Oct,  Nov) 

Dry  Wet  Dry  Wet  Dry  Wet  Dry  Wet 

Deciduous  Forest  2o!o         TOO  33  TOO  22  TOO  4^7  TOO" 
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Coniferous  Forest 

10.0 

70.0 

2.7 

70.0 

2.2 

70.0 

3.3 

70.0 

Wetland/Swamp* 

jO.u 

/U.O 

1  z.  1 

/U.U 

J 1  .J 

/y)A) 

1  Z.v 

/U.U 

Grassland* 

50.0 

70.0 

3.3 

70.0 

3.3 

70.0 

6.6 

70.0 

Cropland* 

50.0 

70.0 

3.3' 

70.0 

4.6 

70.0 

7.9 

70.0 

Urban^ 

10.0 

70.0 

10.0 

70.0 

10.0 

70.0 

10.0 

70.0 

Open  Water 

70.0 

70.0 

70.0 

70.0 

70.0 

70.0 

70.0 

70.0 

Snow/Ice 

50.0 

70.0 

*  in  winter,  wetland,  grassland,  and  cropland  treated  as  a  snow  surface.  bare  soil  and  aciive 

growth. 

"  bare  soil  and  senescent  growth.  ^  consists  of  a  mixture  of  deciduous  forest  and 

buildings. 


Rc  (HNO3):  10  s/m  (for  all  seasons  and  all  surfaces) 

Rc  (HNO2):  10  s/m  (for  all  seasons  and  all  surfaces) 

Rc  (NH3):  28  s/m  (dry  conditions) 

9  s/m  (wet  conditions) 

201  s/m  (when  T2<0°C) 
Rc  (particulates):        0  s/m 

Rc  is  calculated  based  on  surface  wetness  criteria,  such  that  it  either  represents  a  "total 
dry  condition,"  "total  wet  condition,"  or  "weighted  wet  condition"  using  the  following 
flowchart,  and  relative  humidity  (RH)  and  surface  wetness  (SW)  criteria: 
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RH  available 

Go  to  4 
Adapted  from  Bates  (1996) 


Default  =  Rc  value  for  dry  conditions 
Wet  Rc  =  Rc  value  for  wet  conditions 
Weighted  Wet  Rc  =  Time  weighted  wet  Rc 


Time  weighted  wet  Rc 


SW 

(mo 


X  wet  Rc 


SW 

Too 


SW  available 


xdry  Rc 


13 
Go  to  5 


Calculation  of  Vd  in  the  absence  of  meteorological  data: 

Missing  hourly  meteorological  data  are  treated  in  the  following  manner: 

•  1  hour  of  meteorological  data  missing  — >  the  average  resistance  of  the  hours  before 
and  after  are  used  to  represent  the  missing  hour 

•  consecutive  hours  of  meteorological  data  missing  — >  each  hour's  calculated  median 
resistance  for  the  month  is  used  to  represent  the  missing  hours 
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APPENDIX  II 


Environment  Canada  Calculation  Methods  for  Gases  and 

Particulates 

(after  Brooks,  1999a,  199b;  Zhang  et  al.,  2003a,  2003b,  2002a,  2002b,  2001a,  and  2001b) 
Parameters  (for  gases); 


Vd'.  deposition  velocity  (m/s) 

Ra'.  aerodynamic  resistance  (s/m) 

Rb'.  boundary-layer  resistance  (s/m) 

Re'.  surface  resistance  (s/m) 

k:  von  Karman  constant  (0.4) 

u*:  friction  velocity  (m/s) 

u:  horizontal  wind  speed  at  reference  height  (m/s) 

L:  Monin-Obukhov  length 

H:  sensible  heat  flux 

z:  reference  height  ( 1  Om) 

zq:  surface  roughness  length  (m) 

:  stability  function 

Ri:  Bulk  Richardson  number 

g:  gravitational  acceleration  (9.81  m/s^) 

T2:  temperature  at  reference  height  (°K) 

T2p:  potential  temperature  at  reference  height  (°K) 

Tg:  surface  temperature  (°K) 

Td:  temperature  difference  (°K) 

P:  surface  pressure  (kPa) 

Tavg:  temperature  average,  T^^,^  ={T^p+T^)l2 

SR:  solar  radiation  (w/m^) 

RH:  relative  humidity  (%) 

Free:  precipitation  (mm/hour) 

FCLD:  fraction  of  cloud  covering  (%) 

cos  6 :  cosine  of  solar  zenith  angle 

cos  a  :  cosine  of  sun/leaf  angle,  set  as  0.5 

Di:  molecular  diffusivity 

Dv:  water  diffusivity 

Vi:  thermal  diffusivity 

MW:  Molecular  Weight,  MWair  =  29 
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Wsi'  fraction  of  stomatal  blocking  under  wet  conditions 

Rst'.  stomatal  resistance 

Rm:  mesophyll  resistance  (dependent  only  on  chemical  species) 

Rac:  in-canopy  aerodynamic  resistance  (not  chemical  species-dependent) 

Rg:  ground  resistance  (scaled  from  SO2  and  Os's  Rg) 

Rcut:  cuticle  uptake  resistance  (scaled  from  SO2  and  Os's  Rcut) 

Gs(PAR):  unstressed  leaf  stomatal  conductance 

PAR:  Photosynthetic  Active  Radiation 

f{T):  conductance-reducing  effects  of  air  temperature  T  (°C) 

f{D)\  conductance-reducing  effects  of  water-vapour-pressure  deficit  D 

D\  water-vapour-pressure  deficit 

f{y/)'.  conductance-reducing  effects  of  water  stress  y/ 

y/ :  water  stress 

i//^.^  and     2  •  leaf-water-potential  dependency 

E^(T):  saturation  water  vapour  pressure  (kPa)  at  air  temperature  T  (°K) 

E:  ambient  water  vapour  pressure  (kPa) 

f    , :  snow  cover  fraction 

»/  snow 

Sd:  snow  depth  (cm) 

Sdmax-'  maximum  snow  depth  (cm) 

Rcutwo,  reference  values  of  wet  cuticle  resistance,  see  Table  1  for  RcutwO  (O3) 
Pcutwo{^02):     50  s/m  or  100  s/m  for  rain  or  dew  conditions,  respectively 

Rcutdo,  reference  values  of  dry  cuticle  resistance,  see  Table  1 

LAI:  Leaf  Area  Index 

LUC:  Land  Use  Category  (26  LUCs  in  gas  model) 


Deposition  Velocity: 


Aerodynamic  Resistance  (for  surfaces  other  than  water/lake): 

ENVC  adapts  the  Acid  Deposition  and  Oxidant  Model  (ADOM)  formulated  by  Pleim  et. 
al.  (1984)  and  further  investigated  by  Padro  et  al.  (1990). 

R^  =  ^—  [0.74  ln(z  /  Zq)-i//'h],  with  a  lower  limit  of  5  and  upper  limit  of  1000 
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Friction  velocity: 


kii 


In 


zO 


(1  +  4.7/?/) 


ku 


In 


9ARi 


1  +  7.4 


k  In 


Uo 


Unstable,  Ri<0 


Stability  function: 


-4.7- 


0.74x2  In 


l  +  Jl-9 


stable,  0  <  —  <  1 
L 


unstable,  -1<  — <0 
L 


Monin-Obukhov  length: 


L  = 


kHg 


with  a  lower  limit  of  -5  and  upper  limit  of  +5 


Bulk  Richardson  number: 


Ri  = 


Tu 


,  with  a  fixed  value  of  1x10"^^  if  SR>0  and  Ri>0 


Temperature  difference: 

=  T2p  -T^ ,  with  an  upper  limit  of  -10"^°  for  negative  values  and  a  lower  Umit  of  10'^^ 
for  positive  values 
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Sensible  heat  flux: 


uT, 


0.74J  [ln(z/zO)J   [{\  +  4JRiy 
k      I'L  9ARi 


0.74j[ln(z/zO)J   r    l  +  5.3[^/ln(z/zO)]-^ 


Stable,  Ri>0 


Unstable,  Ri<0 


Sub  Layer  Resistance: 


R. 


5  f  Vi 


.  x2/3 


\Di J 

Thermal  diffusivity: 


-4^3/2 


Vi 


145.8  X 10""  r 


+110.4 


Molecular  diffusivity: 


0.00  ir 


Di  = 


where. 


)"'+D"'y 

air  gas  f 


D  =  0.369MW„ +  6.29 


Total  Surface  Resistance: 


1  \-Wst  1  1  ,  1-     •  / 

—  =  1  1  ,  With  a  lower  limit  of  10  s/m 
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Fraction  of  stomatal  blocking  under  wet  conditions: 


0,  SR<200Wm~~ 
(SR  -  200)  /  800       200  <  SR  <  600Wm 
0.5  SR>600Wm-^ 


Stomatal  resistance: 

R^,=  

G^{PAR)f{T)f{D)f{y/)DJD^ 

Water  diffusivity: 


0.00  IT 


7/4 


1  ' 

1 

—  +  

\    air  water  / 

Unstressed  leaf  stomatal  conductance: 

G,,  {PAR)  =  F„„  /  r„  (PAi?,„„ )  +  F^^,  I  r„  (PAR^,^, ) 

F,„=2cos4l-e"°'""™'] 

^shade  ~  ~  ^sun 

r^,{PAR)  =  r^^^^{\  +  bJPAR) 

PAR      -  [^diff^'-^''^'''^  +O.QlR,.^{\A-OALAI)e-'°''     LAI  <  2.5    or    SolarRad  <  lOQwrn'^ 
'''''  ~  1^^/^^^"°'"^'"^  +0.07/?,.,(l.l-0.1LA/V-^°^^    LAI  >  2.5    and    SolarRad  >  lOOwm'' 


PAR 


Rj-^  cosa/ cos 0  +  PAR ^1,^^^     LAI  <  2.5    or    SolarRad  <  lOOwrn'^ 

;,0.8 


\  R'fr  cosa/cosO  + PAR, LAI  >  2.5    and    SolarRad  >  200wm-' 


Conductance-reducing  effects  of  air  temperature  T; 


f(T) 


T-T 


opt 
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T  -T 

where  hr 


T  -T 

opt  rrun 


Conductance-reducing  effects  of  water-vapour-pressure  deficit  D: 

fiD)  =  1  -  b^^jD  ,  with  a  lower  Hmit  of  0. 1  and  a  upper  limit  of  1 

Water-vapour-pressure  deficit: 

D  =  EiT)-E 

Saturation  vapor  pressure: 

17.27(7'-273.16) 


E*(T)  =  0,6me 


7-35.86 


Ambient  water  vapour  pressure: 

E  =  EiT)RH 

Conductance-reducing  effects  of  water  stress: 

where  y/  =  -0.72  -  0.001 35/? ,  with  a  lower  limit  of  0. 1  and  upper  limit  of  1 


Mesophyll  resistance: 

Values  of  Rm  for  all  dry-depositing  species  in  AURAMS  gas-phase  chemical  mechanism 
listed  in  Table  1 . 


In-canopy  aerodynamic  resistance: 


where, 

Raco^  reference  value  for  in-canopy  aerodynamic  resistance  (Table  1) 


Dry  Deposition  IVIonitoring  Method  in  Alberta 


45 


LAI:  Leaf  area  index  (Table  1 ) 


For  some  LUCs,  a  range  of  Raco  values  is  given  to  reflect  the  change  of  canopy  structure 
at  different  times  of  the  growing  season.  Raco  values  for  any  day  of  the  year  based  on 
minimum  and  maximum  LAI  values  given  as: 

.  ,       LAI{t)-  LAIimin)      r  ^        •  i 

Kco  (0  =  Kco  (min)  +  ..  X       (max)  -  R,^.,  (mm)] 

LAI(Max)  -  LA/(mm) 


Ground  resistance: 


I 


a(i)    ^  J3(i) 


where. 


2000  s  m      LUC  1-3  and  Snow  Surface 
200  s  m"'        LUC  4-19,  25  and  26 
500  s  m-'  LUC  20-24 


For  snow  surface,  Rg(03)  adjusted  by  including  a  snow  cover  fraction(  /,„,,„  ): 


1      _  l-2/.„  ^  2/,^ 


R^{0,)     R^(LUC)  R^ 


where. 


sd 


fsnow  =  — — '  note  -  both  /^,,^^^  and  2        have  a  lower  limit  of  0  and  upper 

limit  of  1 

R_  =  2000  s/m 


RASO,)  = 


20  s  m~' 
70(2 -r) 

50 

100 

R^,,  (Table  1) 
mm{2R^,,Ry-''-'-'') 


LUC  1  and  3 
LUC  2,  T  as  Surface  Temperature  (°C) 
LUC  4-26,  Rain 
LUC  4-26,  Dew 
LUC  4-26,  Not  Rain  or  Rain,  T>-\°C 
LUC  4-26,  Not  Rain  or  Rain,  T<-\°C 
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with  a  lower  limit  of  100  s/m  and  upper  limit  of  500  s/m 


For  snow  surface,  Rg(S02)  are  adjusted  by  snow  cover  fraction(  f,„^^): 

^    ^  snow    _|_  f snow 

R^iSOy  R^(LUC) 

where, 

sd 

fsnow  =  — — '  note  -  both  f^„^^  and  2        have  a  lower  limit  of  0  and  upper 

limit  of  1 

/?„-V0(2-r),Tas  °C 


Canopy  cuticle  resistance: 

R 


cutwO 


LAl"'u, 


R. 


IAl"'u. 


mm{2R^.^,,R„,e 


0.2(-l-r) 


Rain  or  Dew 


Not  Rain 


)  T<-VC 


where. 


Rcutwo(03)  listed  in  Table  1 . 

Rcutwo(S02)  =  50  s/m  or  100  s/m  for  rain  or  dew  conditions,  respectively. 
Rcutdo  listed  in  Table  1 . 

with  a  lower  limit  of  100  s/m  and  20  s/m  for  SO2  dry  and  wet  conditions, 
respectively 


For  snow  surface,  Rcut  are  adjusted  by  including  a  snow  cover  fraction  ( ) 
1      l-f  f 

'■    J  snow  _|_  J  snow 


R  R  R 

cut  cut  ^snow 


where 


2000  O3 

70(2 -r)  so. 
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Roughness  length: 


,   .  ,      LAI  (t)- LAI  (mm)      r  ^       ^   .  . 

Zo  (t)  =  Zo  (min)  +  .    X  [z,  (max)  -  z,  (min)J 

LAI{Max)  -  LAI{mm) 


Rain  or  Dew  conditions: 

Condition  = 


{Rain  T  >  213.\5°K  and  FreoOAmm/hr 
[Dew         r>  273.15°/:  and  u\>u\m 


where: 

\.5Coedew 


u 


min 


DQ 

where: 

0.3        FCLD  <  0.25 
Coedew  =  \o2   0.25  <  FCLD  <  0.15 
0.1      0.15  <  FCLD  <  I 

and 

DQ  =  0.0622(1- RH)E  \T),  with  a  lower  limit  of  0.0001 
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Parameters  (for  particulates): 

(Variables  for  computing  Ra  were  not  included) 


deposition  velocity  (m/s) 

aerodynamic  resistance  (s/m) 

Rs: 

surface  resistance  (s/m) 

u*\ 

friction  velocity  (m/s) 

u: 

horizontal  wind  speed  at  reference  height  (m/s) 

8- 

gravitational  acceleration  (9.81  m/s  ) 

T2: 

temperature  at  reference  height  (°K) 

LAI: 

Leaf  Area  Index 

LUC: 

Land  Use  Category  (15  LUCs  m  particle  model) 

SC: 

Season  Category  (5  SCs  in  particle  model,  listed  in  Table  5) 

P : 

density  of  particle 

dp: 

particle  diameter 

C: 

correction  factor  for  small  particles 

M: 

dynamic  viscosity  of  air 

v: 

kinematic  viscosity  of  air 

C: 

Cunningham  slip  correction  factor 

Eb: 

collection  efficiency  from  Brownian  diffusion 

Eim: 

collection  efficiency  from  impaction 

Ein: 

collection  efficiency  from  interception 

Rl: 

correction  factor  representing  fraction  of  particles  that  sticks  to  surface 

So: 

Schmidt  number 

D: 

Brownian  diffusivity 

St: 

Stokes  number 

K: 

Boltzmann  constant  (gm  cm  /s  -°K) 

P: 

surface  pressure  (kPa) 

A: 

mean  free  path  of  air  molecules 

Deposition  Velocity: 

=  y,  +— 

1 

Ra 

+  Rs 

Assumes  S04^"  and  NH/  have  the  same  Vd  ;  and  NO3",  Ca^^,  Mg^^  Na^  and  K^"  have  the 
same  Vd  as  described  in  Cheng  et  al.  (2001). 
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Mass  median  diameter  (MMD)  and  geometric  standard  deviation  (GSD)  for  SO4""  chosen 
as  0.35  |am  and  2.0  pm,  respectively  as  reported  in  Wesely  et  al.  (1985). 

MMD  and  GSD  for  Na^  are  taken  as  5.12  pm  and  2.64  pm,  respectively  as  reported  in 
Ruijgrok  et  al  (1997). 

Gravitational  settling  velocity: 

V  =^ 
18// 

Cunningham  slip  correction  factor: 

C  =  l  +  M(i.257  +  0.4e-"""''") 

Aerodynamic  Resistance: 

Ra  is  computed  with  the  same  approach  as  the  gas  model,  but  the  parameters  are  adapted 
from  Table  4. 

Surface  Resistance: 

=  5  ,  with  a  lower  limit  of  5 

where  £q  is  an  empirical  constant  chosen  as  3  for  all  LUCs 

Collection  efficiency  from  Brownian  diffusion: 

E,=Sc-' 

where  y  lies  between  V2  and  ^/s  with  larger  values  for  rougher  surfaces.  (Table  4) 

Schmidt  number: 

Sc  =  v/D 
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Collection  efficiency  from  impaction: 

'    St  " 


a  +  St 


where, 

a:  constant  varying  with  LUC  (Table  4) 
b:  constant,  chosen  as  2 


Stokes  number: 


smooth  surface,  LUCd>,  9,  12,  13,  and  14 


St 


V 

vegetated  surface,  LUCs  other  than  above 


8 

where  A  is  the  characteristic  radius  varying  with  LUC  and  season  (Table  4) 

Collection  efficiency  from  interception: 

d  ^ 


F    =  — 


,  with  a  upper  limit  of  0.6 


where  A  is  the  characteristic  radius  varying  with  LUC  and  season  (Table  4) 

Correction  factor  representing  fraction  of  particles  that  sticks  to  the  surface: 

/?j  =  exp(-5r'^^) ,  with  a  low  limit  of  0.5 

Mean  free  path  of  air  molecules: 

ju  101.3 


1  =  6.54x10" 


1.818x10"'  P 


v293.15y 


Brownian  diffusivity: 

■indu 
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Dynamic  viscosity  of  air: 

145.8xlO-'r''' 


+110.4 


Kinematic  viscosity  of  air: 

Pair 
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JI 

ii 

JE 

4 
JI 

J- 


7  ?  7  ^  ?  T  ?  T  T  ^  7  ^  ^  7  7  7  ^  "  T  V 
-  V  -  T  V  T  ^  ^  V  -  2   2   -   -   -   -  -  -  o  - 


T  <=  f 


2  2 


8§g8g§§8§S    g    °    g    g  ggggg 

1 1  i  i  1 1 1 1 1 1  1  1  I  I  1 1  i  I  i 

iiiiii§iii  i  8  i  i  iilii 
1 1  i  1 1  i  i  1 1  I  i  i  i  I  i  i  1 1 1 


g  g 

i  I 
i  i 
i  i 


giiii.  !       !  !  !  I  I  !      =  .  =  i 


iiiiii 

IIIIH 


MSI 


I  In  III 


oo^ooooooooooooSooooooooooooooo  o 

I  !l  lis?  i  =° 

stPiiilllyiiiltiiliiiiia. 

I 
I  I 

I 

Jiiitiiiliiliiili 

isigioSiiiiiiiiiiigsiiiiiliililil 


Oh 


E  CO 
E  o 


CO      OO      CO      UD  CO 

in    un    LO    Lo  Lo 


un    uo    Lo    un    Lo  un 


^    ^    ir>    LO  ^ 


^    ^  -^^ 


CNJ      in      CM      LO      uo      CM  CNJ 


CNJ    lo    LO    ;;r:  LO 


CN       LO       LO  LO       LO  LO 


CSJ       LO       CM       LO       LO       CNJ  CNI 


c\j    LO    CM    LO    LO    CNi  eg 


LO      lO      LO  LO 


in    lo    LO    in    CM  CM 

CO      LO      LO      V-      O  o 


in    in    in  CM 


OO     ^      CO     T-  ^ 


OO      ^      CM  — 


in  in  LO  in 
CO    c»    o  ■<- 


T-        T-  ^ 


o  — 


LO      LO      LO      LO      ^  ^ 


CO  ^  CO  T-  ^  ^ 
CD      ^      CD     CD      3,  3, 


^       CD      J=  2f 


C:  Z3 


Q   Q   s   o   o  a 


■<—      CM  CO 


O      •<-     CM  CO 


o 

Q. 
(D 

Q 


APPENDIX  III 


WBEA  Dry  Deposition  Calculation  Methods  for  Gases  (SO2  and  NO2) 

(after  EPCM,  2002,  2000) 

Parameters: 

Vd'.  deposition  velocity  (m/s) 

Ra'.  aerodynamic  resistance  (s/m) 

Rb'.  boundary-layer  resistance  (s/m) 

Re'.  surface  resistance  (s/m) 

0 :  standard  deviation  of  wind  direction  (rad) 

u<  friction  velocity  (m/s) 

u:  horizontal  wind  speed  at  reference  height  (m/s) 

SR:  solar  radiation  (w/m^) 

SW:         surface  wetness 

Deposition  Velocity: 

Vd=  7  ^  r 

{R^+R^+Rj 

Aerodynamic  Resistance: 

4 

— Y     Stable  conditions, night  time  {SR<5w/m^)  or  6<02rad 

U0 
g 

— Y   Unstable  conditions, day  time  {SR  >  5w/ m^)  and  d>0.2rad 


R 


with  an  upper  limit  of  1000 
u  has  a  lower  limit  of  0.36m/s 


Boundary-Layer  Resistance  for  SO2: 

U 

-  with  an  upper  limit  of  1000,  where  = 


U 
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Boundary-Layer  Resistance  for  NO2: 
Surface  Resistance: 

Same  as  Alberta  Environment  method  (Appendix  I),  where  dry/wet  conditions  determined  by 
SW: 

S  W  >0. 1 ,  wet  conditions 
SW  <0.1,  dry  conditions 


Dry  Deposition  Monitoring  Method  in  Alberta 


59 


Dry  Deposition  Monitoring  Method  in  Alberta  60 


m 

« 

m 
m 

m 


• 

•m 


m 
m 

'« 
m 


LIBRARY  AND  ARCHIVES  CANADA 
BIbllotheque  et  Archives  Canada 


3  3286  54121584  0 


